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The problem of simplifying tensor expressions is addressed in two parts. The first part presents 
an algorithm designed to put tensor expressions into a canonical form, taking into account the 
symmetries with respect to index permutations and the renaming of dummy indices. The tensor 
indices are split into classes and a natural place for them is defined. The canonical form is the 
closest configuration to the natural configuration. In the second part, the Grobner basis method is 
used to simplify tensor expressions which obey the linear identities that come from cyclic symmetries 
(or more general tensor identities, including non-linear identities). The algorithm is suitable for 
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implementation in general purpose computer algebra systems. Some timings of an experimental 
^ ; implementation over the Riemann package are shown. 

O 

oo 
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7'. 1 Introduction 

Recently, some attempts to describe algorithms for simplifying tensor expressions have appeared in 
the literature. |l| || |J The tensor symmetry properties and the presence of dummy indices make the 
problem complex. Fulling et al.jlj described an algorithm to enumerate the independent monomi- 
als built from the Riemann tensor and its covariant derivative. They presented explicit tables for 
monomials of order 12 or less in derivatives of the metric. Although explicit bases are presented for 
monomial of order 8 or less, they have not described neither a systematic algorithm to build the basis 
nor a method to simplify generic tensor expressions built using the Riemann tensor. Using the same 
algorithm, Wybourne and Meller^] enumerated the order- 14 invariants. On the other hand, Ilyin 
and Kryukovp| did present an algorithm to simplify general tensor expressions based on the group 
algebra of the permutation group. Though the method is elegant from the mathematical point of 
view, it is inefficient. Quoting the authors: "hardware development is very fast now, and it will be 
possible to solve problems with 11 indices with the help of our program". Dresse presented in the 
second part of his PhD dissertation a new algorithm to simplify tensor expression based on the 
backtrack algorithms for combinatorial objects. |J Most of his effort has been directed to solve the 
"dummy index problem" (described bellow). 

In this work, the problem is addressed in two parts. In the first, the tensor expressions are 
put into a canonical form taking into account the symmetries with respect to index permutations 

1 On leave from Centro Brasileiro de Pesquisas Ffsicas, Rua Dr. Xavier Sigaud 150, Urea, Rio de Janeiro, Brazil. 
CEP 22290-180. Email: portugal@cat.cbpf.br 
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and renaming of dummy indices. In the second part, the Grobner basis method (see Geddes at 
al.|| and refs. therein) is used to address the linear identities that come from cyclic symmetries 
(or more general tensor identities). The algorithm is suitable for implementation in general purpose 
computer algebra systems, since many functionalities of these systems (besides the Grobner basis 
implementation) are required. 

It is known that the names of the dummy indices have no intrinsic meaning, and the presence of 
two or more of them in a tensor expression creates symmetries with respect to renaming. This leads 
to algorithms of complexity 0(n!) where n is the number of dummy indices. Dummy indices difficult 
the definition of tensor rules and the use of pattern matching. One way to solve the dummy index 
problem is to rename them in terms of the index positions and the name of the tensors (see step 9 
of the algorithm). However, the tensor symmetries may change the index positions, invalidating the 
process. This kind of renaming method is invariant if there is a prescribed method to put the indices 
into a canonical position. 

In section 2, an algorithm to put tensor expressions into a canonical form is described. The 
canonical position of the indices is based on definitions 3 to 6. Some parts of the definitions of this 
section involve conventions that can be changed without losing the "canonicalization" character of 
the algorithm. The ideal format for the canonical tensor is 



j> F+ S u A+ B+ C+ Si 



C- B- A- F- 



where F + , Su, A + , B + , C + , Si, C~, B~, A" and F~ represent sequences of indices whose meanings 
are defined in step 7 of the algorithm. The indices of class Si can have contravariant or covariant 
character. They are placed in the midst of the classes that have the character fixed. This ideal 
format is only achieved in special cases, as when the tensor T is totally symmetric or antisymmetric. 
The canonical position of the indices is the one closest to this format, where the notion of "closest" 
is precisely defined. 

Product of equal tensors is addressed by reduction to one tensor of rank equal to the sum of the 
rank of the factors. This reduction process is used in other parts of the algorithm in order to simplify 
the implementation. 

In section 3, some timings of an experimental implementation of the algorithm over the Riemann 
package are presented. Polynomials built from the Riemann tensor are challenges for any algorithm 
to simplify tensor expressions. I believe that special techniques can improve the timings for the kind 
of symmetries of the Riemann tensor, but none has been implemented for the demonstration of this 
section. The implementation in the Maple system|| can easely handle products of three Riemann 
tensors. 

In section 4, the problem of the simplification of tensor expressions under the presence of side 
tensor identities, which come generally from cyclic symmetries, is addressed. The Grobner basis 
implementation of the Maple system is used to accomplish the full simplification. 



2 The algorithm 



Consider the definition of the normal and canonical functions given in Geddes et al. || Let A be a set 
of algebraic expressions which admits a canonical function. Consider the operations of multiplication, 
addition and contraction of tensors as defined in the tensor algebra. Q [ 10| [ 11 ] If a coordinate system 
has been selected, the tensor algebra can be performed through the tensor components. In this work, 
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a tensor expression is any expression written in terms of non-assigned tensor components obeying the 
rules of the tensor algebra whose coefficient factors are members of the set A. Consider Riemannian 
spaces, in which exists a fundamental metric tensor which establish a relation between the covariant 
and contravariant tensor indices. 



2.1 Definitions 

Hypothesis 1: Tensors do not obey any side tensor identities except symmetries with respect to 
index permutations or symmetries with respect to renaming or the inversion of character of dummy 
indices. 

"Symmetries with respect to index permutations" means that the tensors obey one or more 
equations of the kind 

where e = 1 or e = - 1 and 7r(ii • • • i n ) is any permutation of i\ ■ ■ ■ i n . 



Definition 1: Induced symmetry of a sub-set of indices of a tensor. 

The induced symmetries of a sub-set of indices of a tensor are the symmetries that the sub-set 
inherits from the symmetries of the whole set of indices. Pairs of dummy indices are treated as 
independent free indices, hence not permutable. 

For example, the induced symmetry of the first two indices of the Riemann tensor is the skew 
symmetry. The second and third indices have no induced symmetry regardless any contraction 
between the first and fourth indices. 



Definition 2: Equivalent index configurations. 

Two index configurations^ of a tensor T extracted from a tensor product are said to be equivalent 
if one configuration can be put into the other by the use of any of the following properties: 

1. Character inversion of the dummy indices, 

2. Renaming the dummy indices, 

3. Index permutation allowed by the symmetries of the tensors of the tensor product. 



Definition 3: Suppose the tensor T has n indices and let (Ai, • • •, X p ) be a partition of n where 
p is a positive integer less or equal than n. The indices of T can be grouped in disjoint classes 
Ci, ■ ■ C p where a generic class has Aj indices. The indices of each class can be substituted 
with numbers in such way that the indices of class Ci run from (X}" 1 ! Xj) + 1 to £}=i X d . Consider 
all index configurations of the tensor T and let the indices be substituted with the corresponding 
numbers. The configurations are in one-to-one correspondence with the elements of the symmetric 
group 5 n .[0[^3| The following criteria establish an order of decreasing configurations with respect 
to classes C\ to C p : 

a. Smaller value of the position of the first index of class C\ in the tensor T. If the positions are 
equal, consider the position of the next index of class C\. If the positions of all indices of class C\ 
are equal, consider the positions of the indices of the next classes until C p . 

2 The index configuration is the list of indices of the tensor, taking into account the character of each index. 
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b. Smaller value of the first index member of class C\ that appears in the tensor T. If the first 
indices of class C\ that appear in all configurations are equal, consider the next index member of 
class C\. If the indices of class C\ appear in the same order, consider the order of the indices of the 
next classes until C p . 

Definition 3a alone compares the position of the classes, while definition 3b alone compares the 
order of the indices. Given a set of equivalent configurations of a tensor T, definition 3 allows one 
to select the smallest configuration of the set with respect to a given partition of the number of 
indices. The smallest configuration is unique. If definition 3b is not applied, or if it is applied for 
some but not all classes, instead of having one smallest configuration, one may have a sub-set of 
smallest configurations. 

Definition 4: Character normal configurations. 

Let C + and C~ be the classes of the contravariant and covariant indices respectively of a tensor 
T extracted from a tensor product. Consider the set of all equivalent index configurations of T. The 
character normal configurations consist of the sub-set of smallest index configurations according to 
definition 3a with respect to classes C + and C~ . 

Definition 5: Index normal configurations. 

Consider the definition of group I and II given in step 3 and the definition of classes F + , Sn, A + , 
B + , C + , S/|] C~ , B~ , A~ and F~ given in step 7a. The present criteria applies when a tensor (T) 
of group I is extracted from a tensor product. Consider the set of all equivalent index configurations 
of T. Let the indices be relabelled as described in step 7c. The index normal configurations of the 
tensor T consist of the sub-set of smallest configurations that are character normal configurations, 
and satisfy definition 3a for classes F + , Sn, A + , B + , C + , Sj, C~, B~ , A~ and F~ and definition 3b 
for classes F + , S n , B+, C+ Si , C" B~ and F~ . 

Definition 6: Index canonical configuration. 

The index canonical configuration is the only element of the set of index normal configurations 
which fully satisfies definition 3 with respect to classes F + , Sn , A + , B + , C + , Sj , C~ , B~ , A~ 
and F~ in this order. 



The order of the indices of classes A + and A~ and the order of the indices of the sub-classes of Si 
(not including Si °) are not considered in the definition of the index normal configurations. In fact, 
if classes A and Si have less than two elements each, the set of index normal configurations has one 
element, which is the index canonical configuration. 

An algorithm to simplify tensor expressions must recognize when an expression vanishes. The 
following three lemmas guarantee that the algorithm presented here can recognize null tensor prod- 
ucts. 



Lemma 1: Let V be a product of non-null tensors and suppose that no factor vanishes. Let A be 
be a factor and define the tensor S as 

V 

S = T 



3 Si is a collection of sub-classes, which have the same status of the other classes in this definition. The order of 
the sub-classes is described in step 7c. 
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Suppose that A and S share n contracted indices. Let s be the rank of S, and a be n plus the number 
of free indices of A. 

The product V is zero if and only if there exists a factor A having the symmetry 

where the permutation 7r acts only on the indices contracted with the indices of S, and V is invariant 
under the application of the permutation 7r on the corresponding indices of S, that is 

P .V ;::r ..,,, (2) 

where n fs are equal to n i's and the permutation 7r acts on names, not on positions.^ A may have 
indices contracted internally which have not been represented in (1) and (2). Symmetry (1) takes 
account of permutations and character inversions of the dummy indices within A. 
Proof: (=^)(by reductio ad absurdum) There are two cases: 

1. If V is not invariant for any factor A that has the symmetry (1), then 

V^A^S <n ... M . 

Using (1) and renaming the dummy indices, it follows that V ^ —V, therefore 

2. If no factor admits a symmetry ir of the form (1), from the supposition that no factor of V 
vanishes, it follows that V ^ 0. 

So far only products of non-vanishing factors have been considered. What are the conditions that 
cancel a single generic tensor with some or all indices contracted? The answer can be obtained from 
lemma 1. Suppose that T is a tensor of rank m with 2n indices contracted. This tensor can be 
written as 

g g . . T* jl — in jn fl — fm-ln) fO\ 

y lljl y ^njn ' V / 

where a is the permutation of iiji ■ ■ ■ i n jn fi • •• fm-in t na t specifies the actual order of the indices 
of T, and g ij is the metric supposedly symmetric. The free indices are /i • • • f m -2n- For this case, 
the tensor S of lemma 1 is 

^ il jl ' ' ' in jn = 9 hjl 9 injn) 

and is symmetric under the interchange of i p j p into j p i p for p < n, and is totally symmetric under 
the pair interchange of % p j p into i q j q for p,q < n. Tensor S does not vanish (again from lemma 1). 

If the second factor of (3) does not vanish, there are two cases to consider regarding lemma 1. 
In the first case, tensor T has antisymmetry (1) while S is symmetric under the same permutation 
acting on the corresponding indices independently of S being multiplied by T. In the second case, 
the indices of S has the symmetry (2) only if one considers the contractions of indices between S 
and T. These cases reflect items (i) and (ii) respectively of the following lemma. 



4 The character of the indices of A and S need not be contravariant and covariant respectively. The only restriction 
is that the character of the dummy indices are opposite. 
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Lemma 2: Let T be a non-null tensor of rank m with 2n indices contracted (2n < m). Let a be a 
permutation such that 

J 1 jl ••• in jn fl ■■■ fm-2n) 

describes an index configuration with m free indices, such that after j p (p < n) indices are lowered 
by the metric terms as described in (3), one obtains the actual index positions of T. Suppose that 
(4) does not vanish. T is zero if and only if at least one of the following items is satisfied. 

(i) Consider (4). Independently of any contraction, there is a permutation p acting on i\ j\ ■ ■ ■ i n j n 
such that 

rppoa{i\ jl ••• i n j n fl ■■■ f m -2n) _ g p>o(il jl •• • i n jn fl " • fm- 2n ) 

where e = 1 or e = -1 and at least one of the following items is satisfied. 

(1.1) T poa ^ x ^ "' %n Jn ^ '/ m - 2 «) is antisymmetric under one or more interchanges of i p j v into j p i p 
for p < n 

(1.2) T poa ^ %x Jl -^infi ■■ -fm-2n) j s antisymmetric under one or more pair interchange of i p j p into 
i q j g for p,q<n. 

(ii) There is an index character configuration such that T is antisymmetric under a permutation 
7r acting on the contravariant indices (like (1)) and T is invariant under the same permutation 
acting on the corresponding covariant indices (like (2) with all indices in the same tensor). 



Proof: PJ 



For example, suppose that T is a tensor of rank 6 with the symmetries 

rpij klmn rp klij van 

rpi j klmn rp j iklmn 

r-pi j klmn rpijlkmn 

pi j klmn pii j klnm 

Consider the index configuration 

T ikl ikl (5) 

which is equivalent to zero. Due to the contraction of the first and fourth indices, (5) is antisymmetric 
under the interchange of k + and l + , and is symmetric under the interchange of k~ and l~ . This is 
an example of item (ii) of lemma 2. 

There is one case not analyzed yet. A tensor T may vanish due to symmetries regardless of any 
index contraction. This case has no practical applications but the algorithm must recognize what are 
the combinations of symmetries that cancel the tensor. That recognition can be done at the moment 
of defining the tensor, before the algorithm is executed. From now on, suppose that tensors are not 
zero if there is no index contraction. 

Let M. be the set of all mathematically equivalent tensor products generated by all possible 
equivalent configurations of the indices of the tensors of V (V is a product of non-null tensors when 
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all indices are considered free). Definitions 4-6 can be extended to products of tensors by application 
on each factor. 

Lemma 3: Two elements of A4 have the index canonical configuration with opposite signs if and 
only if V is zero. 

Proof: (=>) If V is mathematically equivalent to V and —V simultaneously then V = 0. 

(<=) If V is zero, either one of the factors vanishes or, using lemma 1, there is a factor A 
antisymmetric on some of its indices (like (1)) such that V satisfies (2). If one of the factors vanishes, 
item (i) or item (ii) of lemma 2 are obeyed. For all possible cases, two equivalent opposite terms with 
the index canonical configurations are generated due to the presence of contracted antisymmetric 
indices. 

In worst cases, the algorithm recognizes that a factor is zero in step 7f, that a product is zero in 
step 7h, and that a sum is zero in step 10. 

Rule 1: Consider a tensor product. The character of the first (from left to right) index of a pair of 
summed indices (if some exists) is contravariant, and the character of the second is covariant. The 
same rule applies for the summed indices within a single tensor. 



2.2 The algorithm 

The algorithm is divided into steps grouped by types of action that must be followed in increasing 
order unless otherwise stated. The main goal is to put the indices into a canonical position. The 
canonical position involves the relative position of the indices inside the tensor as well as their 
character. Definition 6 is a precise specification of the canonical position. All dummy indices are 
renamed and the original names are thrown away. The renaming of the dummy indices takes into 
account their position inside the tensors and the order of the tensors. Therefore the renaming process 
leads to canonical names only after the dummy indices have been put into a canonical position and 
an ordering for the tensors has been established. 

Step 1. (Expanding) Expand the tensor expression modulo the coefficients. Select the tensor factors 
of the first term.Q From now on, consider how to put a tensor product into a canonical form. Steps 
1 to 9 are applied to all terms of the expanded expression, one at a time. 

Step 2. (Raising or lowering indices) If there are any metric tensor with indices contracted with other 
tensors, the corresponding indices are raised or lowered according to the contraction character. The 
same goes for other special tensors built from the metric. 

Step 3. (Splitting in group I and II) Split the product in two groups. Group I consists of tensors 
with symmetries. Group II consists of tensors with no symmetries. Tensors of group I are placed in 
the left positions and tensors of group II in the right positions using the commutativity property of 
the product. 

5 For expressions consisting of one tensor, only steps 6 and 9 are performed if the tensor has no symmetries; steps 
7 and 9 are performed if the tensor has symmetries. Classes Sj and Sji are empty is this case. 
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Step 4a. (Merging equal tensors) Tensors of group II with the same name and same number of indices 
and with no free indices merge into a new tensor. Suppose that each original tensor has n indices 
and that there are N equal tensors, then the new tensor has nN indices and is totally symmetric 
under the interchange of the groups of the n indices. After the merging, the resulting tensors are 
incorporated into group I. The information about the relation with original tensors is stored, since it 
will be used at the end of the algorithm to substitute the new tensor with the original tensor names. 

Step 4b. Tensors of group I with same names, same number of indices and same number of free 
indices merge to form a new tensor. The new tensor has the same symmetry under the interchange 
of group of indices as described in step 4a, and each group of indices inherits the symmetries of the 
original tensors. In step 4a, only tensors of group II with no free indices are merged. In this step, 
the tensors may have free indices. 

Step 5a. (Sorting the tensor names) The tensors are lexicographically sorted inside each group. 
Tensors with same name but with different number of indices are sorted according to the number of 
indices. 

Step 5b. Tensors with the same name and same number of indices are sorted according to the 
number of free indices. 

Step 5c. Tensors of group II with the same name, same number of indices, and same number of 
free indices are sorted according to the name of the first free index. Notice that the free indices have 
fixed positions for tensors of group II. 

Step 6. (Fixing the index character of group II) Consider the tensors of group II. The summed indices 
that are contracted within the tensors or those that are contracted with other tensors of group II may 
have their character changed in order to obey rule 1. All other summed indices (the ones contracted 
with tensors of group I) are put covariant. The free indices remain untouched.^. 

Step 7a. (Splitting in classes) Step 7 is performed for all tensors of group I. Consider the first tensor 
of group I (the current tensor). It has symmetries that can involve all indices or only some of them. 
The discussion that follows applies only for the indices involved in the symmetries or contracted with 
indices involved in the symmetries. The pairs of summed indices not involved in the symmetries 
must obey rule 1 and the free indices not involved in the symmetries remain untouched throughout 
the whole algorithm. Consider the indices that can have their positions affected by the symmetry. 
They are split into classes: 

Class F + \ contravariant free indices 

Class Sif. summed indices contracted with tensors of group II. 
Classes A + , B + , C + , A~ , B~ , C~: summed indices inside the tensor 
Class Sf. summed indices contracted with tensors of group I. 
Class F~: covariant free indices 

6 This step can be displaced and performed together with step 9 
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Let / + , sji, a, b + , c + , sj, F + , A, B + , C + , Sj, C~, B~ and F~ respectively. Classes A + and A~ have 
the same size. 

The summed indices inside the tensor are members of classes A + , B + , C + , A~ , B~ or C~ . When 
both indices of a pair of summed indices are involved in the symmetries, the contravariant index is 
a member of class A + , and the covariant is a member of class A~ . If a index of a pair of summed 
indices is involved with the symmetries while the corresponding one is not, there are four cases. Let 
us call i + the index involved in the symmetries and i~ the corresponding index not involved in the 
symmetries. These indices are equal but have different characters (here the signs + and - do not 
describe the character). Suppose that the relative position of i + with respect to i~ cannot be inverted 
due to the symmetries. If i + is at the right of i~, then i + is a member of class B + ; if i + is at the 
left, then i + is a member of class B~ . If the relative positions of i + and i~ can change, then i + is a 
member of class C + or C~ depending on whether i + is contravariant or covariant. 

As soon as class A is determined, it is verified whether there are antisymmetric contracted indices. 
If so, the tensor is null and the algorithm returns to step 2 for the next term.[] 

The indices of class Sj are contracted with the indices of the tensors of group I. Some of the latter 
indices may not be involved in the symmetries. They form sub-class Sj°, which is further split into 
Sj 0+ and Sj° , corresponding to the indices of Sj° contracted with the tensor to the right or to the 
left of the current tensor respectively. The remaining indices of Sj are split into sub-classes. The 
indices of the current tensor contracted with the next (to the right) tensor of group I form the first 
sub-class (Sj x ); the indices contracted with the next tensor form the second sub-class (Sj 2 ), and so 
on until the last tensor of group I. If the current tensor is the first of group I, the sub-division in 
classes is complete; otherwise the sub-division continues and the next sub-class consists of the indices 
the current tensor contracted with the first tensor of group I. The following sub-class consists of the 
indices contracted with the second tensor of group I and so on until the last tensor of group I that 
has not been considered yet. 

Step 7b. (Fixing the index character of group I) In order to obey rule 1, the indices of class Sjj are 
put contravariant. The corresponding indices of tensors of group II have already been put covariant 
in step 6. The indices of class Sj are put contravariant if they are contracted with the tensor at the 
right, or covariant if they are contracted with the tensor at the left of the current tensor. The indices 
of class B + are put contravariant and the corresponding ones are put covariant. The indices of class 
B~ are put covariant and the corresponding ones are put contravariant. The character of the indices 
of classes F~ is maintained throughout the whole algorithm. The character of classes A + , C + , A~ 
and C~ may still change. 

Step 7c. (Ordering and the numbering of the indices) This step establishes the order of the indices 
of classes F + , Sjj, B + , C + , Sj, C~ , B~ and F~ to be used in step 7e. Also, it specifies the numbers 
each index receives in order for definition 3 to be applied. To begin with, let us discuss the order of 
the indices of classes F + and F~ . First, sort the free indices of class F + using their original names. 
The first sorted free index is substituted with the number 1, the second by 2, and so on until the f + 
th index, which is substituted with / + . The same process is performed for the indices of class F~ , 

7 The vanishing of the tensor due to the presence of antisymmetric contracted indices can be a natural consequence 
of the application of steps 7d to 7f together with lemma 3. For the sake of efficiency, it is better to verify the presence 
of antisymmetric contracted indices as soon as class A is determined. 
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which receive the numbers / + .. a~ + 1, • • •, / + .. 

Class Sn consists of summed indices contracted with the tensors of group II. At this point, the 
tensors of group II are ordered. The positions of the indices of these tensors are used as a reference 
to order the indices of class Sn. The first summed index in group II (from left to right) that is a 
member of class Su is the first index of Su. It receives the number f + + 1. The second receives the 
number f + + 2 and so on until / + + Sjj. 

The indices of classes B + , B~ , C + and C~ are contracted with indices that have no symmetries; 
therefore, they have an ordering reference. They follow the same method used for the indices of class 
Sn- 

Class Si can have its sub-classes ordered. The first sub-class is Si° , followed by sub-classes Sj 1 , 
Sj 2 and so on, and the last sub-class is Sj° . The indices of the classes Si 0+ and S/° can be ordered 
since they are contracted with indices not involved in the symmetries. They follow the same method 
used for the indices of class Sn- The indices of the other sub-classes cannot be ordered at this point. 
Also, the indices of class A cannot be ordered at this point. 

Here follows, explicitly, the numbers reserved for each class: 



Class F+ 


1, 


2, / + 






Class Sn 


r 


- + / + + 2, 




f + + su 


Class A+ 


r 


- + 8H + 1,-- 




h .. a+ 


Class B+ 


r 


- .. a+ + 1, •■ 


r 


h .. b + 


Class C + 


r 


" .. 6+ + 1, •• 


r 


" .. c+ 


Class Sj: 


r 


.. C+ + 1, ••• 


f + 


.. si 


Class C- 


r 


' .. s/ + l, •• 


,f + 


.. c~ 


Class B~ 


r 


" c~ + l, •• 


■, r 


- .. b- 


Class A~ 


r 


- .. + •• 


■, r 


.. a~ 


Class F~ 


r 


" .. a~ + 1, • • 




■■ f~ 



Step 7d. (Generating the character configurations) Both indices of a pair from class A can change 
their positions due to the symmetries, but the relative position of some pairs may be fixed. The 
characters of the pairs, that cannot invert the relative position, can be chosen such that they obey 
rule 1. The characters of the remainder indices of class A (let do be the number of pairs of class 
A that can invert their relative position) and the characters of indices of class C cannot be chosen 
a priori. Each pair has two states which are contravariant-covariant and covariant-contravariant, 
corresponding to the characters of the indices. The algorithm generates 2 ( - ao+c ) possible character 
configurations by changing the states of pairs that can invert their relative position. 

If the current tensor is totally symmetric, steps 7d and 7f need not be performed. The canonical 
form can be obtained straightforwardly, avoiding the slowest steps. 

Step 7e. (Applying the symmetries) At this point, more than one equivalent configuration may have 
been generated. The symmetries are applied to all configurations in order to perform the following 
tasks.0 The contravariant indices are pushed to the left positions as far as possible and the covariant 

8 The notation /+ .. /~ means / + + sn + a + + b + + c + + sj + cr + b~ + a~ + f~. The variables a + and a~ are 
equal to a. 

9 The application of the symmetries is a straightforward procedure that can be implemented for each kind of 
symmetry. In the case of the Maple system, tensors can be represented by tables and the symmetries by indexing 
functions, which can perform the tasks described in step 7e. 
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indices to right positions as far as possible. A sign change may be generated if there are antisymmetric 
indices. The state configurations that are members of the set of character normal configurations of 
the current tensor are selected. The symmetries are applied to the selected configurations again in 
order to put classes F + , Sjj, A + , B + , C + , Sj, C~, B~, A~ and F~ as closely as possible in the 
smallest position as prescribed in definition 3a, and after that, the indices of classes F + , Sn, B + , 
C + , Sf°, C~, B~ , and F~ are put as closely as possible in their order as prescribed in definition 3b 
(see step 7c for relabelling). The sub-set of index normal configurations is selected. At this point, 
the character and the positions of all classes have been determined. Only the order of the indices of 
class A and of the indices of the sub-classes of Si (not including Sj°) has not been determined yet. 

Step 7f. (Generating equivalent configurations of class A) In general, the ordering of the indices of 
the sub-classes of Si depends on the ordering of the indices of class A and vice-versa. These classes 
must be ordered together. For now, suppose that all classes Si 1 (i > 0) are empty. In this case, steps 
7a to 7g can be performed for each factor of the product since they are independent of each other. 

The aim of this step is the generation of permutations of class A that preserve the character 
arrangement of the indices. The selected configurations of step 7e are submitted to all possible 
re-orderings of class A (class A + plus A~) allowed by induced symmetry of the indices of class A + 
together with the indices of class A~ . The character normal configurations are selected. 

If the induced symmetry of class A is totally symmetric, the generation of permutation are not 
necessary since the canonical positions can be obtained straightforwardly. 

In many cases, it is sufficient to generate the re-orderings allowed by the induced symmetries of 
class A + and class A~ independently. This kind of re-ordering automatically maintain the character 
configuration. 

Step 7g. (Selecting the index canonical configuration) The indices of all configurations are substituted 
by their correspondent numbers (step 7c), and the index canonical configuration with respect to 
classes F + , Sn, A + , B + , C + , Si, C~ , B~ , A~ and F~ is selected. 

If two elements with the index canonical configuration are selected and they have opposite signs 
then the current product is zero (lemma 3). In this case, the algorithm returns to step 2 for the next 
term. 

Step 7h. (Ordering indices of the sub-classes of Si) The order of the sub-classes of Si and the order 
of the indices of Si have already been determined. The present step finds out the order of the 
indices of sub-classes Si 1 (i > 0) of all tensors of the current product. Consider class A and all sub- 
classes Si 1 (i > 0) of the first tensor. These classes together have a induced symmetry. The same 
can be stated about the other factors of the product. All indices of all classes A and all sub-classes 
Si 1 (i > 0) of all factors are put together in order to form a new tensor totally contracted. The 
order of the indices in the new tensor follows the order of the factors and the order of appearance in 
each factor. The symmetries of the new tensor is composed by all induced symmetries acting in the 
corresponding indices (see example 2). 

By a recursive call of the algorithm, the indices of the new tensor are in class A and are ordered 
by the method described for this class. The dummy indices that come from sub-classes Si 1 cannot 
invert their relative positions (within a pair) hence step 7d need not generate character configurations 
for these indices. 

If the induced symmetry of the indices of the sub-classes Si 1 of a factor coincides with the actual 
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symmetry of these indices (taking into account the contractions of class A) then the indices of class 
A of this factor need not be included in the new tensor. 

Step 8. (Recovering merged tensors) The symmetric (by group of indices) tensors that have been 
formed by merging tensors with equal names, equal number of indices and equal number of free 
indices in step 4a and 4b are converted back by the inverse process to a product of tensors with the 
original names but with the new positions of indices generated by the previous steps. 

Step 9. (Renaming the dummy indices) At this stage of the algorithm, the indices are in their final 
position. The dummy indices are renamed following the rules: There are two cases. The first occurs 
when the whole pair of dummy indices resides inside a tensor. If the tensor name is A and the 
number of indices is m, then the dummy index name will be A_m_i_j, where % is the position of the 
contravariant index and j is the position of the covariant index and _ is some separator .F] If the 
same name appears in other tensors of the product, no name conflict is generated. The second case 
occurs when the dummy indices involve two tensors. Suppose that the first tensor has the name 
A with m indices and the second has the name B with n indices, then the dummy index will be 
renamed A_m_i_B_n-j , where i is the position of the contravariant index inside the tensor A and 
j is the position of the covariant index inside the tensor B. This renaming process proceeds from 
left to right. If a second dummy index receives the same name, then the number 1 is appended in 
its name A_m_i_B_n_j_l , for example. If a third summed index receives the same name, then the 
number 2 is appended in its name: A_m_i_B_n_j_2, and so on.0 

Step 10. (Collecting equal tensor terms) After performing steps 1 to 9 for all terms, collect equal 
tensor products and put the coefficient factors into the canonical form. 

Step 11. (Sorting the addition)^ Each tensor product can be converted to a string by concatenating 
with separators the names of the tensors and the indices in the order they appear in the product. 
These strings are sorted. The order of the terms in the sum is rearranged to be in the same order as 
the sorted concatenated strings. 

2.3 Proof that the algorithm is a canonical function 

Let T be the set of all tensor expressions which obey hypothesis 1. The algorithm described in 
section 2.2 is a function T : T i — > T. 

Theorem: T is a canonical function. 
Demonstration: The proof has three parts. 

1. All operations performed by the algorithm obey the rules of the tensor algebra, hence preserve 
the mathematical equivalence of tensor expressions. 

10 The separator is a symbol not present in tensor expressions. 

11 The method of appending a number to the repeated dummy index names can be fully avoided if step 8 is performed 
after step 9. 

12 In general, this step is not necessary when the algorithm is implemented over multiple purpose computer algebra 
systems. 
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2. For all £1,82 £ T such that £\ = £ 2, ^(£ 1) = ^{£2)- After step 1, £\ and £2 are sums 
of tensor products. Consider a generic tensor product. It is clear from the sorting uniqueness that 
the position of the tensors is unique after the application of the algorithm. The indices of tensors of 
group II and the indices of classes F + , Sn, B + , Si, B~ and F~ of tensors of group I go to a unique 
character configuration, since this is a matter of convention. 

The proof that the indices of tensors of group I (including the merged tensors of step 4) come 
to a unique configuration is consequence of the fact that the algorithm runs over all index character 
configurations of the indices of classes A and C and over all allowed index position configurations 
of classes A and Si. Only one configuration is selected as the canonical configuration unless the 
product is zero (lemma 3). After step 9, the dummy indices have canonical names, completing the 
canonicalization of a tensor product. Steps 10 and 11 put a sum of tensor products into the canonical 
form. 

3. In item 2 is missing the details about the special case when one tensor expression is zero, 
that is, if £ = then T{£) = 0. From lemmas 1 and 2 follow that the indices responsible for the 
cancelation of a factor are in class A and for a product are in class Si. Since the algorithm generates 
the set of all mathematically equivalent tensor products by permuting the indices of classes A and 
Si, lemma 3 garantee that any null product inside the tensor expression is recognized. Other non- 
null tensor products that may still exist are put into a canonical form (item 2) guaranteeing the 
cancellation of a sum of products in step 10. 

In the algorithm, there are shortcuts for special kinds of symmetries avoiding the generation of 
character configurations of step 7d or permutations of step 7f. These shortcuts increase the efficiency 
but they must satisfy items 1-3 of the proof. 



2.4 Examples 

Example 1: Consider the tensor expression R\ a i (contraction of the Riemann tensor). Step 7a 
establishes that classes A and F ~ are the only non empty classes in this case. Class A + is [i + ] , 
class A~ is [i~] and class F ~ is [a, b]. Step 7c establishes that the contravariant index % receives 
the number 1, the covariant index % receives the number 2, and the indices a and b receive the 
numbers 3 and 4 in this order, since a precedes b. Step 7d establishes that there are two character 
configurations to be considered: R l bai and Riba 1 - The symmetries are applied (step 7e) in order 
to put these configurations into the equivalent ones R l bin and R l a ib- Both are selected, since both 
are members of the set of character normal configurations. No further changes are generated by the 
symmetries, so the configurations that are members of the set of index normal configurations must 
be selected. The indices are substituted with their respective numbers, yielding [1,4,2,3] and [1,3,2,4] 
respectively. Definition 3a with respect to the partition (1,1,2) - corresponding to classes [1], [2] and 
[3,4] - selects both configurations. Definition 3b with respect to class [3,4] selects [1,3,2,4] as the 
only member of the set of index normal configurations. Since there are no induced symmetries, no 
more index configurations are generated. The index canonical configuration is [1,3,2,4] . Therefore, 
the canonical form is R % ai b, where i is R_4-l-3, as prescribed in step 9. 

Example 2: Consider the tensor expression T : >S kl T l R i ijk, where S kl is a totally symmetric 
tensor and Rnjk is the Riemann tensor. Step 3 splits the expression into group I: [S kl , Rnj 
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and group II: [T\ T 1 ]. Step 4b merges T J and T l into one totally symmetric tensor (let be called 
T_Ti l ) and adds to group I. Group II is empty now. Step 5a establishes the order for group I: 
[Rnjk, S kl , T_T^ 1 ]. Step 7a establishes that the only non-empty class in this case is Sj. For the 
first tensor one has Sj 1 = [I, k] and Si 2 = [j, i}. The order of these indices has not been established 
yet. These are the only classes since all indices have been covered. Step 7b fixes the character as: 
[R iljk , Ski, T_T ji). Step 7e converts R il i k into R lik i and maintains S k 1 and T_T ji invariant. 
Step 7h generates a new tensor. Let us call JSf 1 tk i ki j i- It has the symmetries of the Riemann 
tensor (excluding the cyclic symmetry) in the first four indices (I + ,i + ,k + ,j + ); is symmetric under 
the interchange of the fifth and sixth indices (k ~,l ~) and is symmetric under the interchange of the 
seventh and eighth indices The algorithm is called recursively, and JV l%k ' ki j i is converted 

to JV 1 1 k i i k i j, determining the order of the indices. After step 8 one has: [R l%k i, Sik, T i, T j}. 
Step 9 establishes that canonical form for the expression is R 1 k J ' S i kT f T j where I = R4 A S-2-1 , 
i = R4_2-T_l_l, k = R4-3S-2-2 and j = R44_T_1_1. 



3 An experimental implementation 

Algorithms to simplify tensor expressions have been implemented in some computer algebra sys- 
tems. [H| flTJ] JET] |L8| Jl9| Some implementations use pattern matching which requires a big database 



of tensor rules, and even worse, sometimes the user must enter the rules. The underlying method 



used by the Ricci package[15] seems similar in some aspects to the method presented here. All these 
implementations have not solved the dummy index problem, therefore the main simplifier spends a 
long time or cannot simplify tensor expressions with many dummy indices. 

In this section, I present an experimental implementation of the algorithm described in section 2 
over the Riemann package. The new package can be obtained from web sites,0 and my purpose is 
to supersede the Riemann package in the near future with the new functionalities introduced to deal 
with tensor components abstractly. 

The function normalform uses the algorithm to put tensor expressions into normal forms. No 
attempt has been made to put the output into the canonical form, that is, ordered with respect to 
the sum and to the product of terms, since this last step is unnecessary for the purpose of any general 
computer algebra system. 

In the examples below, contravariant indices have positive signs and covariant indices have neg- 
ative signs. Tensors are indexed variables and their symmetries are declared by the command de- 
finetensor. Some tensors are pre-defined: Christoffel symbols, Riemann tensor, Ricci tensor and 
Ricci scalar are pre-defined with the names T (Christoffel symbols) and R (Riemann, Ricci and Ricci 
scalar). More details can be found in the new help pages for the commands normalform, definetensor 
and symmetrize, and in the help pages of the Riemann package. 

Expressions that are zero due to tensor symmetries are readily simplified:^ 

> readlib(showtime)(): 

> definetensor(T[i,j,k,l], sym[l,2] and asym[3,4]); 
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See th e addresses: |http : / / www . cbpf . br / ~ port ugal/Riegeom. html or bttp : / / www. astro . qucensu . ca/ ~ port ugal / Ric 



geom.html 



14 All calculations have been performed in a Pentium 120 MHz with 32 Mb of RAM, running Maple V release 5 
over Windows 95. 
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i j k I 



time = 0.01, bytes = 13478 

> exprl := printtensor(R[i,j,k,l]*T[-i,-k,-j,-l]); 

M J- ik j I 

time = 0.05, bytes = 8124 

> normalform(exprl); 



time = 0.12, bytes = 90619 

> expr2 := printtensor(T[i,j,k,l]*V[-i]*V[-j]+V[b]*V[a]*T[-a,-b,l,k]); 

T ijkl V i V j + V b V a T ab lk 

time = 0.01, bytes = 8042 

> normalform(expr2); 



time = 0.30, bytes = 338050 

> off; 

Polynomials constructed with the Riemann tensor exemplify the algorithm's worst performance. 
No special techniques have been implemented for the kind of symmetries of the Riemann tensor. 
In the next example, all possible ways to write the scalars formed by the product of two Riemann 
tensors are generated. The 40.320 expressions reduce to 4 independent non-null forms if the cyclic 
identity of the Riemann tensor is not considered. In the next section, one can verify that the cyclic 
identity reduces to 3 independent scalars (cf. ref. [|ll).[3 

> S := {op(combinat[permute]([a,b,c,d,-a,-b,-c,-d]))}: 

> nops(S); 

40320 

> readlib(showtime)(): 

> SI := map(x->abs(normalform(R[op(1..4,x)]*R[op(5..8,x)])), S); 



{0, 



£> R1R1 R2R2 R3R3 R4R4 ft 



R1R1 R2R2 R3R3 R4R4 



td R1R1 R2R3 R3R2 R4R4 d 

-K it riri R3R2 R2R3 R4R4 



R 



R1R1 R2R2 



R 



R1R1 R2R2 



\{R)\ 2 } 



15 Simplified names for the dummy indices are used for this demonstration. This choice does not provide truly 
canonical names. 
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time = 35094.97, bytes = 21215820858 
The program spends less than one second per expression on average. Next, one example involving 
the product of three Riemann tensors is provided: 

> normalform(R[-c,-d,m,n]*R[a,b,d,c]*R[-n,-m,-b ,-a]); 

T>R1R3 R2R4 R3R1 R4R2 t>R1R3_1 R2R4_1 p 
—tt IX R!R 3 R2R4 JXR3R1 R4R2 R1R3 _1 R2R4-1 

time = 26.41, bytes = 6303158 

> off; 



4 Simplification of tensor expressions 

The algorithm of section 2 achieves a full simplified form of the tensor expressions if the tensors do 
not obey side identities. To accomplish the simplification of tensor expressions obeying side relations, 
the Grobner basis method is used. The general strategy is to put the tensor expression and the side 
relations into a canonical form of the algorithm of section 2, and simplify the new tensor expression 
with respect to the new side relation. 

Here follows an example of how the Grobner basis method is used to simplify the expression 

p a b c d p i I 6) p a b c d p 

> expr := printtensor(R[a,b,c,d]*R[-a,-c,-b,-d]-l/2*R[a,b,c,d]*R[-a,-b,-c,-d]); 

expr:=R abcd R acbd -^R abcd R abcd 

> EXPR := normalform(expr); 

PVpn . p R1R1 R2R3 R3R2 R4R4 p 1 p R1R1 R2R2 R3R3 R4R4 p 

HyArK .— K K RtR! R3R2 R2R3 R4R4 ~2 R1R1 R2R2 R3R3 RW 

> side_rel := printtensor(R[a,b,c,d]*symmetrize(R[-a,-b,-c,-d],cyclic[b,c,d])); 

111 

side_rel := R ab cd (- R ab cd + - R a cdb + - R a db c ) 

> SR := normalform(side_rel); 

nn ._ ^ tdRIRI R2R2 R3R3 R4R4 p 2 j~, RlRl R2R3 R3R2 R4R4 p 

<->-tt ,— 3 R1R1 R2R2 R3R3 R4R4 ~ ^ R1R1 R3R2 R2R3 R ^ RJ > 

> simplify(EXPR,{SR=0}); 



From this example one notices that the scalars R abcd R acbd and R abcd R abcd are not inde- 
pendent. 
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5 Conclusion 



An algorithm to simplify tensor expressions based on computable definitions have been described. 
It has two parts. In the first part, the expression is put into a canonical form, taking into account 
symmetries with respect to index permutations and the renaming of dummy indices. The definition 
of the canonical form involves some conventions that can be changed without disqualifying the 
definition. The conventions are based on the implementation simplicity and on a readable display 
for tensor expression. In the second part, cyclic identities or more general kinds of tensor identities 
are addressed through the Grobner basis method. The expression and the side relations are both 
put into a canonical form for the Grobner method to work successfully. 

In this work, a precise definition of the canonical form for tensor expressions is provided. No 
restriction is imposed on the kind of symmetries that the tensors can obey. For most of the symmetries 
that occur in practise, the algorithm is very fast. The symmetries of the Riemann tensor reveal the 
algorithm's worst performance, but even in in this case it is useful for practical applications. 

The invariant renaming of the dummy indices plays an important role in the efficiency of the 
algorithm, since it neutralizes the symmetries that come from dummy index renaming. This is a 
solution for the dummy indices problem mentioned in the introduction. 

An experimental implementation of the algorithm over the Riemann package 0] is available free 
from web sites.0 All calculations and timings presented in this work can be reproduced in the Maple 
system. 

Acknowledgements : 

I thank Ray McLenaghan and Keith Geddes for pointing out references and for the kind hospitality 
at the University of Waterloo. This work was made possible by a fellowship from CAPES, Ministerio 
da Educagao e do Desporto, Brazil. 



References 

[1] S. A. Fulling, R. C. King, B. G. Wybourne and C. J. Cummins, Normal forms for tensor 
polynomials: I. The Riemann tensor, Class. Quantum Grav. 9 (1992) 1151-1197. 

[2] V. A. Ilyin and A. P. Kryukov, ATENSOR - REDUCE program for tensor simplification, Com- 
puter Physics Communications 96 (1996) 36-52. 

[3] A. Dresse, PhD thesis, Universite Libre de Bruxelles, 1993. 

[4] B. G. Wybourne and J. Meller, Enumeration of the order- 14 invariants formed from the Riemann 
tensor, J. Phys. A25 (1992) 5999-6003. 

[5] G. Butler and C. W. H. Lam, A general backtrack algorithm for the isomorphism problem of 
combinatorial objects, J. Symb. Comput. 1 (1985) 363-381. 

16 See footnote 12. 



17 



K. O. Geddes, S. R. Czapor and G. Labahn, Algorithms for Computer Algebra, Kluwer Academic 
Publisher, 1992. 

Portugal, R. and Sautii , S., Applications of Maple to General Relativity, Computer Physics 
Communications 105 (1997) 233-253. 

Waterloo Maple, Inc., 450 Phillip Street, Waterloo, Ontario, N2L 5J2, Canada. See the web 
pages: |http: / / www.maplesott.com"7| 

D. Lovelock and H. Rund, Tensors, Differential Forms and Variational Principles, John Wiley 
& Sons, 1975. 

L. P. Eisenhart, Riemannian Geometry, Princeton University Press, 1966. 

T. Y. Thomas, Tensor Analysis and Differential Geometry, Academic Press, 1965. 

D. E. Littlewood, The Theory of Group Characters and Matrix Representations of Groups, 
Oxford University Press, 1950. 

B. L. van der Waerden, Algebra, Frederick Ungar Publishing Co., 1970. 
N. Pelavas and R. Portugal, to appear. 

J. M. Lee, D. Lear and J. Roth, Ricci - A Mathematica package for doing tensor cal- 
culations in differential geometry (User's Manual version 1.2) 1995. See the web pages 
[http: / /www.math.washington.edu/~lee/Ricci . 

L. Parker and S. M. Christensen, MathTensor, A System for Doing Tensor Analysis by Com- 
puter, Addison- Wesley, 1994. 

R. Bogen et al., Macsyma Reference Manual, vol II, chapter V2-3, Laboratory for Computer 
Science, MIT, 1983. 

L. Hornfeldt, STENSOR user guide and STENSOR reference manual, University of Stockholm, 
Institute of Theoretical Physics, 1988. (See also M. A. H. MacCallum and J. Skea, Sheep: A 
computer algebra system for General Relativity, Lecture Notes from the First Brazilian School 
on Computer Algebra, vol. 2, eds. M. J. Rebougas and W. L. Roque, Claredon Press, Oxford, 
1994.) 

[19] M. Kavian, R. G. McLenaghan and K. O. Geddes, MapleTensor: Progress Report on a New 
System for Performing Indicial and Component Tensor Calculation using Symbolic Computa- 
tion, ISSAC'96, Proceedings of the 1996 International Symposium on Symbolic and Algebraic 
Computation (ETH 1996), Y. N. Lakshman, Ed., pp.204. 



18 



